%T1
%计算金地距离
clc;clear
format long
N=1000;%插值点数

%原始数据
x0=18:2:30;  %日期
y0=[9.9617724,9.9543645,9.9468069,9.9390950,9.9312245,9.9231915,9.9149925];

x1=18:(30-18)/N:30;
%x1=18:30;
y1=lagrangeinterpolation(x0,y0,x1);
y_exact=9.935799;
%找出离y_exact最近的点
min=1000;
for k=1:N+1
    if abs(y1(k)-y_exact)<min
       min=abs(y1(k) -y_exact);
       m=k;%离y_exact最近的点的序号
    end
end

day=floor(x1(m));%天
hour=floor((x1(m)-day)*24);%小时
min=floor(((x1(m)-day)*24-hour)*60);%分钟
s=floor((((x1(m)-day)*24-hour)*60-min)*60);%秒

fprintf('时间为：%d 天 %d 时 %d 分 %d 秒',day,hour,min,s)
error=abs(y1(m)-y_exact)











